Geometrical Expression for the Angular Resolution of a Network of 

Gravitational- Wave Detectors 



Linqing WerQ 

International Center for Radio Astronomy Research, School of Physics, 
University of Western Australia, 35 Stirling Hwy, Crawley, WA 6009, Australia 

Yanbei Cheid 

Division of Physics, Mathematics, and Astronomy, Caltech, Pasadena, CA 91125, USA 

We report for the first time general geometrical expressions for the angular resolution of an 
arbitrary network of interferometric gravitational-wave (GW) detectors when the arrival-time of a 
GW is unknown. We show explicitly elements that decide the angular resolution of a GW detector 
network. In particular, we show the dependence of the angular resolution on areas formed by 
projections of pairs of detectors and how they are weighted by sensitivities of individual detectors. 
Numerical simulations are used to demonstrate the capabilities of the current GW detector network. 
We confirm that the angular resolution is poor along the plane formed by current LIGO- Virgo 
detectors. A factor of a few to more than ten fold improvement of the angular resolution can be 
achieved if the proposed new GW detectors LCGT or AIGO are added to the network. We also 
discuss the implications of our results for the design of a GW detector network, optimal localization 
methods for a given network, and electromagnetic follow-up observations. 

PACS numbers: 04.80.Nn, 95.75.-z, 97.80.-d, 97.60.Gb 



I. INTRODUCTION 

Several types of astrophysical sources are expected to 
be detectable both in gravitational waves (GWs) and in 
conventional electromagnetic (EM) wavelengths. For ex- 
ample, long gamma-ray bursts have been conjectured to 
originate from asymmetric core collapse of massive stars, 
and short gamma-ray bursts might be produced by the 
coalescence of compact binary objects containing neutron 
stars. Both of these could emit gravitational waves in 
the frequency band of ground-based laser interferometer 
GW detectors (e.g., Ref. [T]). Several large-scale inter- 
ferometric GW detectors have reached (or approached) 
their design sensitivity, and are coordinating to operate 
as a global array. These include the LIGO detectors 
at Livingston, Louisiana, and Hanford, Washington, US, 
the Virgo detector in Pisa, Italy, the GEO 600 detector 
in Hannover, Germany, and the TAMA 300 detector in 
Tokyo Japan. Upgrades to existing detectors (Advanced 
LIGO and Advanced Virgo) have been planned [lHl]j 
while new detectors (LCGT in Japan [F and AIGO in 
Australia [6]) are still being proposed. In case of a strong 
EM event, follow-up searches for GW signals can be con- 
ducted in archived data in the time window of the event 
(e.g., Ref. [7]). On the other hand, EM follow-ups to 
probable GW events require a clear understanding of the 
angular resolution of an array of GW detectors. 

The angular resolution of an individual GW detector, 
arising from its antenna beam pattern, is rather poor [S] . 
However, the large baselines of the current GW-detector 



network facilitate better angular resolution via triangula- 
tion. Several localization algorithms have been proposed 
and the effect of arrival timing uncertainties as well as 
amplitude information of GWs have been investigated [S]- 
[Hj. Quantitative studies of the angular resolution of a 
network of GW detectors have been conducted by several 
authors, both for a ground-based detector network and 
for the future space GW detector LISA pDH25] . A stan- 
dard approach is to calculate numerically the Fisher in- 
formation matrix, which leads to a method-independent 
lower bound on the statistical errors of estimated param- 
eters (see a review in Ref. [2T). On the other hand, 
explicit analytical expressions for the network angular 
resolution are rare in the literature largely because of 
the complexity involved in derivations. 

Two approximate analytical expressions for the angu- 
lar resolution can be found in the literature (summarized 
in [25]) for a network of three GW detectors. One is an 
elegant approximate geometrical formula for 3 detectors 
due to Thorne (as cited in Eq.(8.3) of Ref. 9J): the solid 
angle uncertainty is 



An = 



2e^Ari2Ari3 
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where c is the speed of light, Ati2 and Ati3, are time- 
of-arrival accuracy between pairs of detectors, A is the 
area formed by the three detectors, and 9 is the angle be- 
tween the source direction and the normal to the plane 
of the three detectors. However, the underlying assump- 
tions and derivation of this expression are not available 
in the literature. The dependence of angular resolution 
on the signal-to-noise ratio (SNR) was derived by Tinto 
in Ref. ^£,, by expressing the above time-of-arrival accu- 
racy as a function of SNR and frequency [26] (derived 
from the Fisher matrix assuming all other information of 
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the waveform as perfectly known). The other formula is 
based on the numerical result of the angular resolution 
of the LIGO- Virgo network around a wave incident di- 
rection normal to the plane formed by the three detector 
sites P3 — for GWs emitted from neutron star-neutron 
star (NS-NS) inspirals [22] • This particular resolution 
was then rescaled by the cosine of the wave incidence an- 
gle and SNR [25,. Analytical geometrical expressions or 
approximate ones for the angular resolution for an array 
of more than three detectors have not been obtained in 
the literature. 

In this paper, we deduce explicit analytical expressions 
for the angular resolution of an arbitrary GW detector 
network in terms of observables such as cross-sectional 
areas of the network and energy flux of the incoming 
GW. We use only the time-of-arrival information, ignor- 
ing additional (usually rather poor) information from the 
directional-derivatives of antenna beam pattern functions 
— and therefore arrive at a conservative estimate. Such 
an approximation allows us to obtain expressions that 
have explicit geometrical meanings, further generalizing 
Thome's formula to an arbitrary number of detectors, 
and several particular scenarios. In particular, we con- 
sider both short signals, during which motion of the de- 
tector network is negligible, and long signals, for which 
the trajectory traced by the detectors during the signal 
determines the effective size of the detector network. We 
also consider signals with known or unknown waveforms, 
but always assume unknown arrival times of the signal. 
In this paper, the scenarios where signals have completely 
known and completely unknown waveforms are termed 
interchangeably as the best-case and the worst-case sce- 
narios respectively. 

We focus on deriving explicit expressions for several 
situations that will arise in the practice of searching 
for and localizing GWs. Specifically, we derive general 
expressions for short signals assuming only that arrival 
time is one of the unknown parameters (summarized in 
Eq. (41) and text thereafter) and for long signals assum- 
ing known waveform (Eq. (49 1). Based on these formulae. 



cuss the implications for the design of a detector network 



we show simplified solutions for Eq. (41) in several real- 
istic situations: (1) short signals in the worst-case sce- 



nario (Eq. (26l) and in the best- case scenario (Eq. (42)) 
for an arbitrary network of detectors, (2) special cases 
of the two- and three- detector networks in the best case 
(Eq. (31 1-([32|) and their representations when the wave 
is short and monochromatic (Eq. (33)-(34])), (3) long sig- 
nals in short and long observations with detectors in cir- 



cular motion (Eq. (51) and Eq. (52) respectively). 

This paper is organized as follows. In section |lTj we 
explain our notation. In section |III[ we derive analyti- 
cal expressions of the angular resolution for an arbitrary 
detector network and for special cases. We show explicit 
derivations for the worst-case and the best-case scenario 
in section IIII A II and section IIII A 21 and then derive a 
general expression in section [III A3| In section |IIIB[ wc 
derive a general expression for long-duration wave and its 
application to detectors at circular motion. We then dis- 



and localization strategies in section IV The astrophys- 
ical applications of our results are shown in section [V] 
In section |Vlj we discuss the possible errors in our esti- 
mation caused by ignoring the directional derivatives of 
antenna beam patterns. Our results are summarized in 
section IVIIl 



II. MATHEMATICAL PRELIMINARIES 

A. Antenna Pattern of a Detector Network 

In this paper, we assume a network of Nd gravitational- 
wave detectors, with spatial locations given by the vec- 
tor r/ with / = 1,2,... each with spatial size much 
smaller than the GW wavelength. The strain of an in- 
coming GW observed by an individual detector / is then 
a linear combination of the two wave polarizations in the 
transverse traceless gauge, 



di[to + Ti+t) 



f+h+{t) + f^h^{t), Q<t<T, (2) 



where Iq is the arrival time of the wave at the coordinate 
origin and r/ is the time required for the wave to travel 
from the origin to reach the J-th detector at time t, 



Ti{t) = n-vi{t)/c. 



(3) 



Here n is the propagation direction of a GW, t a [0, T] 
is the time label of the wave, and T is the signal du- 
ration. The quantities /+ and are the detector's 
antenna beam pattern functions |27j for the two wave 
polarizations (/i+, h^). They depend on the relative 
orientation between the detector configuration and the 
frame in which the polarizations are defined (which is in 
turn related to the source direction — n). In particular, 
given a Michelson-type interferometer with orthogonal 
arms along e^, and ej,, and given the symmetric, trace- 
free polarization tensors e+(n) and ex(n) to which the 
wave polarizations refer, we have 

/+,x (n) = (e^ (g) - (g) By) ! G+^x (n) (4) 

where the symbol : stands for contraction. Note that dif- 
ferent conventions can be used to define g+^x (n), as long 
as these tensors are symmetric, trace-free, and satisfy 



n- G+,x(n) = 



(5) 



The Fourier transform of the time-series data from the 
/-th GW detector is 



di{n) = / di{t)e'"*dt . 



(6) 



Denoting the corresponding one-sided noise spectral den- 
sity by Si(fl), we define a whitened data set in the fre- 
quency domain. 



di{n) = Sj '-{n)di{n). 



(7) 
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Vector d(rj) then corresponds to the whitened data set 
at each frequency. For short-duration signals where mo- 
tion of the detector array is unimportant, antenna beam 
patterns are treated as constant, hence 



Ah{n) , 



(8) 



where ^ is a Nj^ x Nj_ diagonal matrix with 0/j = Sjjflrj, 
or 



cf) = n 



(9) 



A is an Nd x 2 matrix of the antenna pattern functions 
for all detectors weighted by noise. 



fi (") 



(10) 



and h{Q) is a 2-dimensional vector function containing 
the two polarizations of a GW in the frequency domain. 



(11) 



For simplicity, we keep the ri-dependence in the nota- 
tion only when it is necessary for clarity. At any indi- 
vidual frequency fl, even though the response vector d 
is Nd dimensional, the existence of only two independent 
signal polarizations, -I- and x , means that the set of all 
possible signal vectors is 2-dimensional. 



B. Fisher Matrix 

We define the following inner product between two vec- 
tors. 



(a|b> 



2 / — a^^b 



2tt 

+ 00 



2^ 



a*j{n)bi{n). 



(12) 



Under the assumption of stationary Gaussian detector 
noise, the optimal squared signal-to-noise ratio is j28j . 

P% = (d|d). (13) 

The Fisher matrix for parameters 9 can be defined as 

r,, = {de,d\ds,d) ; (14) 

note that data are already whitened. The Cramer-Rao 
bound states that for an unbiased estimator (the en- 
semble average of which is the true value), the Fisher 



matrix sets a method-independent lower bound for the 
covariance matrix of estimated parameters when consid- 
ering statistical errors. In the case of high signal-to-noise 
ratio, the covariance matrix of a set of parameters 6 is 
approximately given by the inverse of the Fisher matrix 
(see [m [23l [30l [31] and references therein). 



V.J = {T-%. 



(15) 



Now suppose that in addition to 0, there exist more un- 
known parameters, A. Then the Fisher matrix pertaining 
to 9 can be written as 



min(W|5d) = -F.jC 



,69, + 0{\69\^) , (16) 



where 



Jd = d( 



,X + SX) -d{9,X) . 



(17) 



Using this formulation, A can also consist of a contin- 
uum of parameters — the minimization would become 
a quadratic variational problem. Note that for a dis- 
crete parameter set, the covariance matrix can also be 
calculated using standard matrix inversion for a matrix 
consisting of block matrices. In this paper, we show our 
derivations based on the variational method. 



III. ANALYTICAL EXPRESSIONS OF 
ANGULAR RESOLUTION 



We calculate the angular resolution of a detector net- 
work by applying the Fisher matrix to obtain method- 
independent lower limits [29 on the statistical errors in 
estimating the direction of a GW source. The limits are 
for unbiased estimators and for Gaussian noise (for cau- 
tions in using these limits, see |30j). 

We denote Z as the direction of n, and then the er- 
ror in solid angle (measured in steradians) is defined by 
covariance of nx and ny as 

A^s = 271 ^J{Snj^){Snl.) - {SnxSny}^ = 27r(det T)-^/^ , 

(18) 

where 6nx and Sny are the deviations of nx and ny 
from their true values and T is the Fisher matrix per- 
taining to angular parameters nx and ny. The factor of 
2tt is introduced so that the probability that estimated 
parameters fall outside an area Ai7 due to statistical error 
is e^^^/^^= . For angular parameters in a polar coordi- 
nate system with colatitude 9 and longitude (/>, we have 
Afl, = 27r|cos6i|^(A6'2)(A(/>2) - {A9A(j))^. In the next 
two subsections, we derive expressions for det T and Aflg 
for short and long signals. 
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Short-duration GWs 



where 



1. Worst-case scenario: signal with unknown waveform 

For a short signal with unknown waveform, the un- 
known parameters consist of the sky coordinates of the 
gravitational- wave source Uj ( j — 1, 2) and the unknown 
waveform h(rj). Varying all unknown parameters, we 
have (note that Sto is absorbed into Sh), 



<.^djct)kh{n)5nj +A(5h(f2) 



(19) 

where we have ignored the change of A induced by Srij 
(see sec. VI for a discussion of its effect) and defined 



(9,0) ^ 



dcj) 

dn-i 



(20) 



In the following, recall that d, (5d, h, (5h, (p and A are 
frequency dependent, but we drop their frequency depen- 
dence in equations for simplicity. We then have 



((5d|5d) = ^{{dj4>)kh.\{dk(f))A.h)5nj5nk 

3 

- iY,{^5h.\{dk<j>)A.\i)6nk 

k 

+ (A(5h|A(5h). 
Minimizing ((5d|(5d) over Jh, we obtain 

6\i ^iY^{k^k)-^k\dkcj))k\i5nk , 



(21) 



(22) 



and 



r,fe = ((a,0)Ah|(i-p)(afe0)Ah) 
= ((a,0)d|(i-p)(9fc0)d) 



(23) 



where we have defined I as the identity matrix and 

P = A(A^A)-iA^ , (24) 

which is the projection operator into the whitened sig- 
nal space, i.e., the 2-dimensional space spanned by the 
columns of A. 

After a straightforward algebraic manipulation, we ob- 
tain 



[detr] 



8c4 



jif AiMKr^j X tml) ■ np (25) 



AjK = {ndj\iSjK - P.JK)ndK) 

= iSjK - Pjk) f 2 j —O'd^dA , (26) 

and rjK = - r,/. 

Note that |(rxjxrML)-n|, is just twice the area formed 
by the projections of the detectors J, L and M onto 
the plane orthogonal to the wave propagation direction. 
The quantity /S.jk can be interpreted as the projection of 
the weighted data correlation between detectors J and K 
into a null space. This means that when the waveform is 
unknown, source localization is possible only when there 
exists a null space, which in turn must remain null when 
we get the propagation-direction right. This property can 
be used for source localization [S] and for a consistency 
check for a GW from a given direction [TUl [TH [52] . In ad- 
dition, the value of Pjk is directly related to the inverse 
of A^^A. Thus if A is singular or ill-poscd in the sense 
that the one of the singular values is much smaller than 
the other, a pseudo-inverse should be used to calculate 
the correct Fisher-matrix. In terms of implementation, 
this corresponds to ignoring data corresponding to weak 
network sensitivity [TH [151 133] 



2. Best-case scenario: signal with known waveform hut 
unknown arrival time 



Now let us turn to the best-case scenario with known 
waveform (but unknown arrival time), we have 



5d = -le-'^ 



Ah. 



(27) 



A similar derivation leads to 

r,fe = {{djcf>)d\{dk4>)d) 

_ ((9,-(/>)d|r!d)(f]d|(9fc(/.)d) 
(r2d|rjd) 

After some algebraic manipulation, we obtain 

_ J_ Cj6f^L6/|(rifj X vml) ■ np 

8c4 ^ 



(28) 



J,K.L,M 



1^1=1 



(29) 



Here we have defined 



2 / ^n'\dj?. 
, Ztt 



(30) 



J,K,L,M 



Note that O is directly related to noise-weighted energy 
flux received by the detector (c.f. Ref. [8]). 

Two-detector configuration in the best-case scenario: 
In this case, only a 1-dimensional angular parameter can 
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be determined. The 1-cr 1-dimensional angular resolution 
of a 2-detector network can therefore be derived directly 
from the trace of the Fisher matrix in Eq. (28), with 



error angular separation defined to be twice the standard 
deviation AO^ = 2/^tr{T), 



^^(2, best) _ 



2c 



(31) 



where D±^ is the projected distance of the two detectors 
in the plane perpendicular to the wave direction, and the 
best determined direction is normal to the equal time- 
delay lines of the two detectors. The factor of 2 accounts 
for the total width of the l-cr angular separation. We 
leave the scaling of Ci,^2r ■ ■ by their sum to empha- 
size the importance of the fractional contributions of the 
noise-weighted energy flux from individual detectors. 

Three- detector configuration in the best-case scenario: 
The simplified expression for a 3-detector network de- 
rived from Eqs. (flSl), (pOl) and pOl) is 



^Q(3,best) 



TTC 



1 



(6 + 6 + V ^1^26/(^1 + 6 + 6)3 ' 

(32) 

where A± — |(ri2 x ris) • n\/2 is the projected area of 
the network perpendicular to the wave-propagation di- 
rection. Again, the fractional contributions of the noise- 
weighted energy flux coupled to individual detectors play 
an important role in the angular resolution. 

For the special case of a monochromatic GW at fre- 
quency /, for a 2-detector network. 



best) 



PN fDj_ TT Y pIpIIp% ' 



For a 3-detector network. 



best) 



1 



pI PA^ AttV pIpIpI/p%' 



(33) 



(34) 



where pj is the optimal SNR-squarcd for detector / and 
pj^ is the optimal network SNR-squared [Eq. (13l]. 

Compared with previous expressions summarized in 
Sylvestre [ 



discrete set of unknown parameters represented in a vec- 
tor A of size Nx. Using the same variational approach 
described previously, or a standard matrix inversion tech- 
nique, we obtain the Fisher matrix for the angular reso- 
lution 

r,, = (a,d|a,d)-(a.d|%d)(a,~d|a,^d)"'(%d|9,d) 

= (S.dll-Plcl.d) (35) 

where P = \d^d){d^d\dxdf\dxd\, d^d is the Na x Nx 
derivative matrix of the data d with respect to the un- 
known parameters A with the matrix components de- 
fined as {d^d)ii — ddi/dxi and the bracket opera- 
tion of two matrices is defined as a matrix (A|B)ij = 

Note that P is a projection operator on the vector 
space spanned by d^d and that I — P is a projection 
operator onto the space orthogonal to that of d^d. Evi- 
dently, we have Pja^jd) = |%d) and (I - P)\d^d) = 0. 

The matrix component of the Fisher matrix can be 
written from Eq. (|35|) as. 



(36) 



Ty = y^^AjjdiTjdjTj, 



where 

Aij = {ndj\ndj)Sij - {ndi\dxJi)B-l{dxAj\^dj) 

where djj is a Kronecker delta 



1 for I = J 
for / 7^ J 



(37) 
(38) 



and B = {d^d\d^d) is a matrix with components 
B„ = 2l ^dx^d^dx.d. 



(39) 



For a general assumption that the arrival time to is 
one of the unknown parameters, A/j = 0. To 
see this, assume that to is the fc'-th unknown param- 



Eq. (34) includes a new feature indicating eter. Then Aj, 



to, dx„di 



-Old I for any / and 



that for a given network SNR, the angular resolution is 
limited by the least-sensitive detector. When one detec- 
tor has null response to the wave (e.g., pi ~ 0), we have 
AQ 3> 1, until limited by contributions from the direc- 
tional derivatives of the antenna beam pattern functions 
we ignore in Eq. ( 19 1. This is expected because if the net- 



Ik' — ~i^j{c>\idj\0.dj) for any column number I. 



work contained only 2 detectors then only one dimension 
of the source direction could be resolved. 



3. Short signals: general expression 

Previous derivations can be extended to a general case, 
where in addition to angular parameters, there exists a 



.1 

= Y,{^di\ndj)dij - Y,{ndi\dxJi)B]:l Y.{dxjj\^dj) 

J kl J 

= (ndiindi) -iJ2{^di\dxJi)B-i^Bik' 

kl 

- {ndj I ndj) - J2 {^di I dx, dj)6k 



>kk' 



{fldi\ndi) -i{ndi\dx^,di) 
0, 
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where we have used ^)kiBik 
enables us to obtain from Eq. ( 36 ) 



-'fcfe 



This proof 



[dot r] 

J,K,L.,M 



AjA'ALAfKrxj X tml) ■ np, (40) 



In summary, for an arbitrary network of GW detec- 
tors and arbitrary incoming GWs, the component of the 
Fisher matrix for angular parameters is given by Eq. ( 36 1 . 
With the general assumption that wave arrival time to is 
one of the unknown parameters, the angular resolution 
has a simple geometrical form, 



the same expression as Eq. (25) in sec. IIIAl 



^^(Short) ^ 



4^/27^c2 



T,j,K,LM ^jk^lmHykj X Vml) ' nj^ 



(41) 



where A/j is defined in Eq. (37) for the general case, in write 



Eq. (26) for the worst-case scenario, and it can be shown 
that for the best-case scenario, 



A 



ij 



for Ij^J, 



(42) 



consistent with previous results in Eq. ( 29 ) and Eq. ( 30 ) 



Note that only J ^ J terms contribute to the angular 
resolution. That is, once the waveform is known, the 
angular resolution does not depend on the correlation of 
the data between detectors. 



B. Long-duration signal with known waveform 

Finally, when a long signal with known waveform con- 
tinues for a long time, it is the area mapped out by the 
trajectory of the detector network that determines the 
angular resolution. In order to see this, we first note 
that there exist two time scales in our problem: motion 
of the detector network (slow), and the signal and detec- 
tor response time scale (fast). We first rewrite 6d from 
Eq. ^ 



5di{t) 



-5to ^Y.'^\{t)5T 



t- to- 



ri • Vi{t) 



-6to-J2 



r\{t) 



di{t) 



(43) 



Here we have assumed that the antenna pattern changes 
at a frequency much lower than that of the signal, and 
that the speed of the network is much lower than the 
speed of light. As before, we have also ignored the change 
of antenna beam patterns induced by 5n. We can then 



{5di\5di) = / Ut) 
Jo 



-St.-Y^'^Sn, 

A — ^ r 



dt (44) 



where 

iiit) = 2 J di{t - T/2)wi{T)di{t + T/2)dT , (45) 

where wi is the inverse Fourier transform of 1/5'/ ($7). 
Note that ^/ can be viewed as the rate of increase of 
the signal-to-noise ratio at time t — when we regard d 
as signal, with the original noise spectrum. In obtain- 
ing Eq. (44), we have assumed the observation time, as 



well as the time scale at which r/(<) changes, to be much 
longer than the signal correlation time (i.e., the range 
of integration for r). Following the same variational ap- 
proach described previously, the Fisher matrix can be 



obtained by minimizing Eq. (44) over (Stoj we obtain 



1 r 

-2 V,rk 



J 



(46) 



where ^ is defined in Eq. ((30^ 



Tk 



(47) 
(48) 



are the average of the fc-th coordinate of detectors in the 
network, throughout detection time, and correlations be- 
tween the j-th and k-th coordinates. When the detector 
trajectory has a size much bigger than that of the array 
(as in the case of observing GWs from pulsars using the 
existing LIGO- Virgo network), we can omit the J index 
of Tj, simply replacing it by the mean position of the 
detector network r*^, and the weighted average is simply 
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over time. The determinant of the inverse of the Fisher 
matrix yields the angular resolution, 



EjO(r)Y/var(re)Var(r^) - Cov\re,r^) 

(49) 

where 9, (j) represent the longitude and latitude like coor- 
dinates (e.g., right ascension and declination in the celes- 
tial coordinate system). For a general polar coordinate 
system, rg = dgn ■ r and that = d^n ■ r. Var and Gov 
are variance and co- variance with expectation values cal- 
culated according to Eqs. (47) and (48 1. The angular 



resolution is now clearly related to the square of the area 
mapped out by the trajectory of the network. 

Detectors on circular orbits: As a simple example, we 
consider the situation where the detector network makes 
a circular motion with radius i?* and angular frequency 
w*. Supposing that are all constants, then in the situ- 
ation where the network trajectory has mapped out a size 
much larger than the size of the network, yet w*T 1, 
we find 



[Var(re)Var(r'0) - Cov^(re,r^)] = 



8640 



i?f sin i"?, , 



(50) 

where i„ is the angle between the source direction and 
plane of the circular motion. The factor sin^ i„ accounts 
for the projected area from by the detector network's 
trajectory. Therefore, for a short observation. 



^^(Co„t,S) 



12V15c- 



(51) 



where pT is the optimal signal-to-noise ratio within the 
observation time T. In this case, the error area decreases 
like T-'^ (c.f. Refs. [26ll34|), where T'^ comes from the 
increment of SNRs with time and (w*T)~^ from the in- 
crement of the area formed by the trajectory of the entire 
detector network. For longer observations of u;*T ^ 1, 
contributions from the area saturate because the maxi- 
mum area is ^ 7ri??, 



^^(ContX) 



Pp'^TTRllsini 



(52) 



That is, the error area decreases like T^^ as a result of the 
increment of total SNR-squared. Multiple detectors can 
be treated in the same way as a single detector, since the 
area formed among several detectors is much less than the 
area formed by the detector's trajectory over the typical 
observation time of much longer than minutes. Multiple 
detectors thus contribute mainly by increasing SNR. 



IV. IMPLICATIONS 



We emphasize the clear geometrical meaning of | (r^ j x 
Tml) • n| in Eq. (41 1 for short signals, which is twice the 



area formed by the projections of detectors J, K, L and 
M onto the plane orthogonal to the wave propagation 
direction. This means that, when all detectors from the 
network form a plane, the angular resolution is better for 
sources at directions perpendicular to the plane but is 
poor for directions along the plane. A similar geometri- 
cal term is shown in Eq. ( 49 ) for long signals of known 
waveforms. 

Our formula is consistent with the understanding that 
a larger network improves angular resolution. As we show 
in sec. V, addition of either LCGT or AIGO to the cur- 
rent ground-based LIGO- Virgo detector network results 
in a significant improvement of the angular resolution for 
short signals. Inclusion of the southern hemisphere de- 
tector AIGO will improve dramatically the angular res- 
olution of the network and also break the degeneracy in 
the angular resolution along the plane formed by detec- 
tors on the northern hemisphere [6l I35j . For observations 
of long signals where the size of the network is less im- 
portant than the area mapped out by the trajectory of 
the network, long observing time is essential. 

For short signals, besides geometrical area, the angu- 
lar resolution of a detector network can be improved by 
maximizing values of A/j in Eq. (41 ). In the worst-case 



scenario where the waveform is unknown, the angular res- 
olution largely depends on the noise- weighted correlation 
of responses between detectors (Eq. (26l). In the best- 
case scenario, when wave parameters except for arrival 
time are known, the best achievable angular resolution 
strongly depends on the fractional energy flux coupled 
to individual detectors and does not depend on correla- 
tion of data between detectors. Therefore, building de- 
tectors of parallel antennae and comparable sensitivity 
is advantageous for localizing GWs in either improving 
the data correlation (worst case) or balancing fractional 
energy fluxes among detectors [best case). 

In the best-case scenario, the angular resolution is ap- 
proximately inversely proportional to the total received 
GW energy flux weighted by noise. When the noise spec- 
tral density Si{f) can be replaced with a characteristic 
value, Si{f) ^ So, it follows from Eq. (32) that for a 
3-detector network, 



oc 



So 



AE, 



GW 



(53) 



where AEgw = ^1 Ef=i 1^ /o" f^//'M/(/)l' is the 
GW energy coupled by the network (c.f. Ref. [S]). For 
GWs from coalescing binaries of two neutron stars (in- 
spiral source), even though the mergers occur around 1.5 
kHz, the noise- weighted energy flux peaks around 150 
Hz as this is where the LIGO detectors are the most 
sensitive. This makes the inspiral GW sources relatively 
low-frequency sources in terms of localization. 

For a given detector network of fixed configuration, 
the same considerations should be given for localization 
strategies. For sources of known waveforms, we need to 
focus on maximization of the total coupled energy flux 
weighted by noise while maintaining a balanced budget 
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for fractional energy flux coupled to individual detectors. 
For localizing GWs of unknown waveforms, it is crucial 
to maximize the noise- weighted data correlation between 
detectors and to make use of the maximum dimension of 
the null-space O [TH [H [33] . 



are only two detectors available at different sites, say the 
two LIGO detectors (L and H), we can locate the source 
with 1-dimension angular resolution using Eq. (33). The 
angular separation of the 1-cr error-bars is 



V. ASTROPHYSICAL APPLICATIONS 

We first apply our formulae to the existing detector 
network, that is, the 3-detector network consisting of the 
LIGO detectors at Livingston (L) and Hanford (H) in the 
US and Virgo (V) in Italy (the network is often called 
the LHV network) with 4 km-armlengths. When there 



^^(Sliort,2f) 



, 150Hz 10 3000km 

f Jn D ' 



I (1/2) 

pIpI/p 



N 



/2/2 



sm Ir, 



(54) 

where D is the distance between the two detectors, i„ 
is the angle between the 2-detector baseline and wave 
direction. The error bar is about a factor of 3 better for 
the LIGO- Virgo baseline. For a 3-detector network, 



^^(Short,3f) ^ 8 sq_degs 



150Hz 10 

/ PN 



r 



1/27 V2/2 



An pIpIpIIPn I 



(55) 



where Ajq is the triangular area formed by the three 
detector sites in the detector network. We adopted 
Am = lO^^cm^ for the LIGO- Virgo network. i„ is the 
angle between the wave direction and the plane formed 
by the three detectors. 

We also apply Eq. (511 to the ground-based observa- 
tions of monochromatic GWs, e.g., from pulsars. For 



a typical observation time longer than minutes, the an- 
gular resolution is determined by the area of the vir- 
tual network formed by the detectors' Earth-Sun motion 
which dominates over the area formed by Earth's self- 
rotation. The "&est case" scenario for observations of 
minutes <C T ^ months is 



^^(Cont.S) _ 2.4 sq-arcmins 



/900Hz 10 
\ f PT 



V2/2 



(56) 



For longer observations, our approximation of monochro- 
matic wave is most likely no longer valid as higher-order 
derivatives of frequencies generally must be considered 
(e.g. Ref. [26]). 

The angular resolution for a monochromatic source 
with detectors on the Earth-Sun orbit can also be writ- 
ten using Eq. (52 1 for long observations, e.g., more than 



a year's observation using the future space detector LISA 
at its low frequency range. 



^^(Cont.L) _ Q 73 sq_degs 



/3mHz30y 72/2 
\ f PT 



smir, 



(57) 



Our result is consistent with a simple estimation using 
the difi^raction limit [35] but is a factor of a few smaller 
than the results in Ref. [H] probably because our result 
represents the best-case scenario where we assume per- 
fectly known waveforms. 

We show in Figs. [T]|7]sky maps and statistical behav- 
ior of the error ellipse of angular parameters derived for 



short signals in Eqs (41 1 for the existing three detector 



network of LIGO and Virgo and for a possible 4-detector 
network in the future. With the assumption of station- 
ary Gaussian noise, at 95% confidence level, the area of 
the error ellipse 



(58) 



Fig. [T] shows an all-sky map of the error ellipses in the 
best-case scenario calculated from Eq. (28) for the LIGO- 



Virgo 3-detector network at the advanced configuration 
to detect model GWs from the bar-mode instability of 
supernova (the Ull model as it is from Ref. 03)- The 
duration of the wave is about 20 ms with central fre- 
quency around 600 Hz. Similar results are obtained if we 
use more sophisticated waveforms from Ref. [55]. The 
GWs are injected uniformly in 270 sky directions at a 
given time. The distances to the source are adjusted so 
that the 3-detector optimal network signal-to-noise ratio 
Pn — 10. Fig. [2] and Fig. [3] show the same all-sky maps 
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FIG. 1: All-sky map of error ellipses of angu lar parameters 
at the 95% confidence calculated from Eq. ( 28 1 for the 3- 
detector network of the Advanced LIGO detectors (L and 
H) and the Advanced Virgo (V) to detect GWs from bar- 
instability that center around 600 Ifz for the best-case scenario 
[see also Eqs. (41 1,(42 l]. These ellipses were calculated at a 
fixed time for a fixed optimal network SNR pM ~ 10. Shown 
also in the background are contours of light arrival-time delays 
between detector pairs at a 2 ms interval for the L-H pair and 
4 ms intervals for all other pairs. 




FIG. 3: Same as Fig. [T] but with AIGO (A) from Australia 
added to the LIGO- Virgo network. For AIGO, the Advanced- 
LIGO sensitivity is assumed. 
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FIG. 4: Cumulative distribution of areas of error ellipses for 
angular parameters at 95% confidence, Af2o.95, for the best- 
case (shown in Fig. [l] and in Fig.[2| and corresponding worst- 
case scenario for the LIGO- Virgo (LHV) and LIGO-Virgo- 
LCGT (LHVC) network. 



FIG. 2: Same as Fig. [T] but with LCGT (C) in Japan added 
to the LIGO- Virgo network. For LCGT, the Advanced-LIGO 
sensitivity is assumed. 



if the planned GW detector LCGT (C) or AIGO (A) 
is added. The noise spectral density of all detectors are 
drawn from the design sensitivity for the Advanced LIGO 
detectors [55^ (with zero-detuning of the signal recycling 
mirror and high laser power). 

It is apparent that for the best- case scenario, the shapes 
of error ellipses are determined by contours of equal light 
arrival-time delays between detectors. With addition of a 
new detector to the 3-detector network, there is a signifi- 
cant improvement to the angular resolution due to longer 
baselines and improved SNRs. Fig. [T] shows that when 
all detectors are in the northern hemisphere, the angu- 
lar resolution is relatively poor along the plane formed 



by these detectors. Compared to LCGT, the location of 
AIGO in the southern hemisphere more efficiently lifts 
the plane-degeneracy formed by detectors in the north- 
ern hemisphere. 

Fig. |4] and Fig. [5] show cumulative distributions of 
Ar2o.95, the areas of error ellipses shown in Figs. [l]-[3] 
The improvement with a larger network at the worst- 
case scenario is more prominent than the best-case one 
because, in addition to improvement in SNRs and ge- 
ometrical area, the dimension of the null-space doubles 
when the fourth detector is added (sec III A 1 ) . Our re- 



sults show that for the advanced configuration, 50% of 
the sources with SNR pjy = 10 can be best localized 
within 8 sq-degs for the best cases and within 50 sq-degs 
for the worst cases for the LIGO- Virgo network. These 
areas will be reduced by a factor of 2.5 and 10 for the best- 
case and the worst-case scenarios respectively if LCGT 
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FIG. 5: Same as Fig. [4| but with AIGO added as the 4th 
detector (denoted the LHVA network) instead of LCGT. 
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FIG. 6: Cumulative distribution of the areas Afio.gs of error 
ellipses at 95% confidence for the heat- and for the worst- 
case scenario for GWs from coalescing NS-NS binaries at the 
optimal network SNR pat = 15 at random times and random 
orientations of the binaries (see text). We use the 3-detector 
network of the initial LIGO detectors (L and H) and initial 
Virgo (V) at the design sensitivity. 



FIG. 7: Cumulative distribution of fractional difference of 
the angular resolution for GWs from coalescing binaries of 
neutron stars with and without contributions from direction 
derivatives of antenna beam pattern functions. Negative val- 
ues in X-axis correspond to overestimation of Af^o.gs. Same 
data for Fig. |6] are used. 
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FIG. 8: Comparison of cumulative distributions of Af2o.95 for 
the LIGO- Virgo network for GWs from the bar-mode instabil- 
ity with and without contributions from direction derivatives 
of antenna beam pattern functions. 



is added to the LIGO- Virgo network (Fig. [4]) . Similarly, 
factors of 6 and 15 reduction are observed for these er- 
ror areas for the best-case and the worst-case scenarios 
respectively if AIGO is added instead (Fig. [s]) . 

For the interesting case of GWs from coalescing bina- 
ries of two neutron stars, we show in Fig.|6]the cumulative 
distribution of the solid angle AJIq.qs of the error ellipse 
for the network consisting of the initial LIGO (LH) detec- 
tors and initial Virgo with the design sensitivity [121 lU] • 
We adopt the conventional second-order post-Newtonian 
approximation for the phase of GW waveforms and New- 
tonian amplitude [121 US]- The masses of the neutron 
stars are chosen to be 1.4 Mq each. The distances to 
the source are again adjusted so that the optimal net- 
work SNR is pn = 15. The GW signals are injected 
randomly in 500 sky directions at randomly chosen times 



within a sidereal day and with polarization angles as well 
as binary inclinations drawn from uniform distributions. 

It is apparent in Fig. [6]that 50% of the sources detected 
by this initial LIGO- Virgo network have Ar2o.95 ^ 23 sq- 
degs {best case) and Afio.gs < 80 sq-degs for the worst 
case. This and results for higher frequency burst-like 
GWs discussed in Fig. [4] and Fig. [5] mean that wide-field 
EM cameras with reasonable sensitivity and angular reso- 
lution are advantageous for follow-up observations to best 
catch the electromagnetic counterparts of GW sources 
with current LIGO- Virgo network. We have also found 
similar cumulative distributions for the advanced detec- 
tors for the same fixed SNRs. Within such large error 
ellipses, if improved localization (e.g., by X-ray counter- 
parts) is not available, careful strategies need to be de- 
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FIG. 9: Cumulative distribution of fractional difference in 
the angular resolution for GWs of the bar-mode instability 
caused by ignoring direction derivatives of antenna beam pat- 
tern functions. The same data as in Fig. [8] are used. 



signed to identify host galaxies using wide-field optical 
or radio telescopes [3S]. A detailed study in this aspect 
using the most up-to-date galaxy catalog and applying 
the cataloged B-band luminosity to trace possible GW 
events is on-going 44J. The addition of LCGT or AIGO 
can bring a factor of > 10 reduction for the medium val- 
ues of Af2o.95 in the worst-case scenario and a factor of a 
few for the best-case scenario, and will help significantly 
the pointing of EM telescopes for follow-up observations 
and elimination of confusion sources. 



VI. CONTRIBUTIONS OF ANTENNA BEAM 
PATTERNS 

Our derivation of geometrical expressions for the angu- 
lar resolution of a network of GW detectors does not take 
into account the directional derivatives of antenna beam 
pattern functions (Eq. (|l9|)). It is foreseeable that such 
derivatives can be important for low-frequency signals 
if |09er| < \dg\og f^'^l is approximately true, where r 
is time-delay between pairs of detectors. On the other 
hand, this additional effect will mainly lower the lower- 
limits on AQ calculated in this paper, meaning, most of 
the Afig or Afio.ga given here are over-estimates of their 
true values. We define the fraction changes in Afio.gs as 



5(Ar!o.95)/Ar!o.95 - 1 



^^^0.95/ A^a95' 



(59) 



where Afig.gs is the error area we calculate in this paper, 
AOj^gg is the actual value if the derivatives of antenna 
beam pattern functions with respect to sky directions are 



taken into account in Eq. ( 19 ) 



The inspiral sources, one of the most important GW 
sources for the ground-based detectors, are considered to 
be "low-frequency" sources as discussed previously with 
the noise- weighted energy flux peaked around 150 Hz for 
a 1.4-1.4 Mq binary GW source for the initial or Ad- 
vanced LIGO. The effect of the directional derivatives are ignored (see sec. VI I. In reality, careful designs, e.g.. 



of antenna beam pattern functions is therefore not com- 
pletely negligible compared to the time-delay effect. We 
found that (Fig. [6]) the cumulative distributions with and 
without such effect are consistent with each other within 
22% for Ar2o.95 l£ 23 sq-degs in the best-case scenario. 
The discrepancy is significantly larger (> 60% ) for the 
worst-case scenario where Ario.95 > 35. That is, the 
effect on the cumulative distribution of Afio.QS is signifi- 
cantly larger for large values of Ar2o.95 and for the worst 
cases. On the other hand, Fig. [7] shows that, in the best- 
case scenario, about 35% of Ar^o.gs are overestimated by 
more than 50%. For the worst-case scenario, about 62% 
of Ai7o.95 are overestimated by more than 50%. Fig. [t] 
also shows that a much smaller fraction of AJlo.gs are 
underestimated. 

The effect of the directional-derivatives of antenna 
beam pattern functions for GWs at higher frequencies 
is demonstrated for the supernova bar-instability case in 
Fig. [8] and Fig. |9] In this case, the central frequency 
is about 600 Hz. As a result of higher frequencies, com- 
pared to the inspiral sources, the difference in cumulative 
distributions of Aflo.gs is less significant (Fig. |8]). The 
difference is less than a few percent for the best-case sce- 
nario and less than 10% for the worst case scenario for 
Ar2o.95 < 100 sq-degs. Fig. |9] shows that 15% of Al^o.gs 
are overestimated by more than 30% for the best-case sce- 
nario, while 25% are overestimated by the same amount 
for the worst-case scenario. The fraction of Af2o 95 un- 
derestimated are higher than the inspiral sources though. 

While directional-derivatives of the antenna beam pat- 
tern functions contribute to changes of the areas of error 
ellipse, they do not significantly alter the direction of the 
axes. A histogram using the same data for the inspiral 
case shows that 99.6% of the time, the changes of the 
angles are within ±10%. For the example of the burst- 
like case, 98% of the time, the fractional changes of the 
angles of axes are less than ±5%. 



VII. CONCLUSIONS 

In summary, we have derived for the first time ana- 
lytical expressions for the angular resolution for an arbi- 
trary network of gravitational-wave detectors where the 
directional derivatives of antenna-beam pattern functions 
are ignored. Our results demonstrate the explicit depen- 
dence of the angular resolution on the geometrical config- 
uration of the network, the total noise-weighted energy 
fiux coupled to the network, its fractional distribution 
to individual detectors, and correlation of data between 
detectors. These dependences are intrinsic to the con- 
figuration of the detector network and the sources. The 
results are method-independent and they correspond to 
the best possible localization any unbiased methods can 
achieve in the presence of statistical errors. Our estimate 
is conservative especially for low-frequency sources as the 
directional derivatives of antenna beam pattern functions 
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to remove multiple local minima or maxima caused by in- 
terference when combining data from different detectors, 
or to break the mirror degeneracy in arrival time delays 
for the three-detector case (by using the wave amplitude 
information) , are required for localization methods in or- 
der to achieve the "best" limits. 

Derivation of the angular resolution for short signals 
including those for a network of two and three detec- 
tors are presented in Eq. (31 1 and Eq. (32). Our results 



are consistent with what was previously known from the 
diffraction limit: that a larger network yields better an- 
gular resolution. We confirm that the angular resolution 
is poor along the plane formed by current LIGO- Virgo 
detectors and is better for directions normal to the plane. 
Numerical results are included to show how a new detec- 
tor in Japan (LCGT) or in Australia (AIGO) can dra- 
matically improve the angular resolution of the existing 
network (Fig. |4] and Fig. [s] and discussions thereafter), by 
contributing to longer baselines, additional energy flux, 
extended null signal space, and by breaking the plane- 
degeneracy formed by three detectors. 

Compared with previous approximate expressions for 
3-detectors (sec. [l]), our results are more rigorous and 
show the explicit roles of individual detectors. The an- 
gular resolution of a detector network depends on the 
fractional energy flux coupled to individual detectors in 
the best case and on correlation of data between detectors 
in the worst case. In the best-case scenario, the angular 
resolution is largely limited by the least sensitive detec- 
tor. These have significant implications for the design of 
a detector network and for the design of an optimal local- 



ization method for a fixed network (see sec. IV for more 



discussion on implications). Moreover, our results apply 
to an arbitrary network of any number of detectors. 

We have also derived a geometrical expression for long 
signals for a simplified case where waveforms are known 
(Eq. (49)). The situation where detectors are in circu- 
lar motion and the signal is monochromatic is discussed. 
The dependence of the angular resolution on the latitude 
of the source is apparent in our formulae. We also demon- 
strate that the angular resolution improves rapidly with 
the observing time ~ T'* initially with short observations 
and saturates to ~ T for longer observations (Eq. (|5ll) 



Eq. (52)), consistent with previous knowledge [26, .34j7 



We have further presented through numerical simu- 
lations the distribution of the areas of error ellipses at 
95% confidence level for two of the most important GW 
sources for ground-based detectors. (1) We apply our 
calculations to GWs from coalescing binaries of neutron 
stars using the sensitivity curves of detectors that are 
operating at this writing. The actual limit of the angu- 
lar resolution for these inspiral sources should be closer 
to the best-case scenario since theoretical waveforms are 
known and essential parameters can be estimated with 
great accuracy independent of source direction determi- 
nation |45J. (2) We apply our method to burst-like GWs 
using a representing waveform from bar-instability of 
neutron stars in supernovae for advanced detectors. The 
actual angular resolution for this type of "burst" source 
fits in our worst-case scenario for short signals. We shows 
that, for the existing LIGO- Virgo detector network, as- 
suming uniform distribution of sources, at an optimal 
network SNR of around 15, 50% of inspiral sources can 
be located within 23 sq-degs {best case) at the 95% confi- 
dence level. For the burst source, without any knowledge 
of the waveform, at SNR of 10, 50% of the sources can 
be localized within 50 sq-degs (worst-case), but it can 
be reduced to 8 sq-degs if we have predicted waveforms 
available (e.g., from Ref. [37l[38]). Results for the ini- 
tial or advanced detectors are similar. Our results imply 
that, for prompt follow-up electromagnetic observations 
directly using triggers from current GW network, wide- 
field telescopes are desirable. 
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